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THE FOKAS METHOD TO THE SASA-SATSUMA 
EQUATION ON THE HALF-LINE 

JIAN XU AND ENGUI FAN* 



^^ . Abstract. We present a Riemann-Hilbert problem formalism for 

, . the initial-boundary value problem for the Sasa-Satsuma(SS) equa- 

Qh. tion: on the half-line. And we also analysis the global relation in 



this paper. 



1. Introduction 



Several of the most important PDEs in mathematics and physics 
Ch I are integrable. Integrable PDEs can be analyzed by means of the In- 

verse Scattering Transform (1ST) formalism. Until the 1990s the 1ST 

K*" ■ methodology was pursued almost entirely for pure initial value prob- 

^O ■ 
00 ' lems. However, in many laboratory and field situations, the wave mo- 

t:;^- I tion is initiated by what corresponds to the imposition of boundary 

'sj" . conditions rather than initial conditions. This naturally leads to the 

fr^ ■ formulation of an initial-boundary value (IBV) problem instead of a 

pure initial value problem. 

^^ , In 1997, Fokas announced a new unified approach for the analysis of 

^ . IBV problems for linear and nonlinear integrable PDEs [H |2] (see also 

[3]). The Fokas method provides a generalization of the 1ST formalism 

from initial value to IBV problems, and over the last fifteen years, 

this method has been used to analyze boundary value problems for 

several of the most important integrable equations with 2x2 Lax 

pairs, such as the Kortewegde Vries, the nonlinear Schro dinger, the 

sine-Gordon, and the stationary axisymmetric Einstein equations, see 
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2 J.XU AND E.FAN 

e.g. [U |9]. Just like the 1ST on the hne, the unified method yields an 
expression for the solution of an IB V problem in terms of the solution of 
a Riemann-Hilbert problem. In particular, the asymptotic behavior of 
the solution can be analyzed in an effective way by using this Riemann- 
Hilbert problem and by employing the nonlinear version of the steepest 
descent method introduced by Deift and Zhou [15]. 

It is well known that the nonlinear Schrodinger(NLS) equation 

iqT + ^qxx + \q\'^q = o (i-i) 

describes slowly varying wave envelopes in dispersive media and arises 
in various physical systems such as water waves, plasma physics, solid- 
state physics and nonlinear optics. One of the most successful among 
them is the description of optical solitons in fibers. But, by the ad- 
vancement of experomenal accuracy, several phenomena which can not 
be explained by equation (11. ip have been observed. In order to under- 
stand such phenomena, Kodama and Hasegawa proposed a higer-order 
nonlinear Schrodinger equation 

iqr + ^qxx + \q\'^q + ieiPiqxxx + /?2|gP?x + /^s^d^Hx} = 0. (1.2) 

In general, equation (11.21) may not be completely integrable. How- 
ever, if some restrictions are imposed on the real parameters /3i, /32 and 
/Ss, then we can apply the 1ST to solve its initial value problems. Until 
now, the following four cases besides the NLS equation itself are konwn 
to be solvable: 

• the derivative NLS equation-type I(/3i : /?2 : /33=0:1:1), 

• the derivative NLS equation-type II(/3i : (^2 : /33=0:1:0), 

• the Hirota equation(/3i : /32 : /33=1:6:0), 

• the Sasa-Satsuma equation(/3i : (^2 : /33=1:6:3). 

iqr + -qxx + \q\'^q + i^iqxxx + Q\q\'^qx + ^qi\qf)x) = o (1.3) 

Recently, Lenells develop a methodology for analyzing IBV problems 
for integrable evolution equations with Lax pairs involving 3x3 ma- 
trices pJJ . He also used this method to analyze the Degasperis-Procesi 
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equation in p^. In this paper we analyze the initial-boundary value 
problem of the Sasa-Satsuma equation on the half-line by using this 
method. The 1ST formalism for the initial value problem of the Sasa- 
Satsuma equation has been obtained in |10j . 

According to [lOj we introduce variable transformations, 

u{x, t) = q{X, T) exp{^(X - -^)}, (1.4a) 

be 18£ 

t = T, (1.4b) 

X = X - ^. (1.4c) 

Then equation (II. 2p is reduce to a complex modified KdV-type equation 

ut + e{u^xx + Q\u\'^Ux + ?>u{\u\^)x} = 0. (1.5) 

Organization of the paper. In section 2 we perform the spec- 
tral analysis of the associated Lax pair. And we formulate the main 
Riemann-Hilbert problem in section 3. We also analysis the global 
relation in section 4. 

2. Spectral analysis 
The Lax pair of equation (11.51) is [lOj, 



^x = U^, ^ 






(2.1a) 



where 



and 



^t = 1/^. (2.1b) 



U = -ikA + Vi. (2.2) 



V = -Aiek^A + V2 (2.3) 
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In the following, we let e = 1 for the convenient of the analysis. 

2.1. The closed one-form. Suppose that ^(x, t) is sufficiently smooth 
function of (x, t) in the half-line domain Q = {0 < x < oo,0 < t < T} 
which decay as x — ?■ oo. Introducing a new eigenfunction /i(x, t, k) by 



"^ = fie 
then we find the Lax pair equations 



iAkx—4iAk^t 



//^ + [ikA, fj] = Vifi, 
fit + [4zPA, ft] = V2fi. 



the equations in flA.2p can be written in differential form as 

where W{x, t, k) is the closed one-form defined by 

W = e(''="+''^'*)^(yirfx + V2dt)fi. 



(2.6) 



(2.7) 



(2.8) 



(2.9) 
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2.2. The /ij's. We define three eigenfunctions {/Uj}i of (1A.2I) by the 
Volterra integral equations 

fi^ix,t,k)=I+ f e^-'''''-^'''''^^Wj{x',t',k). J = 1,2,3. (2.10) 

J-yj 

where Wj is given by (12. 9p with fi replaced with /ij, and the contours 
{lj}i ^^6 showed in Figure 1. The first, second and third column of 



T' 






T 



>ix,t) 



.{x,t) 



{x,t) 



O 



O 



O 



71 



72 



73 



Figure 1. The three contours 71,72 and 73 in the (x, t)— domain. 



the matrix equation (I2.10p involves the exponentials 






2ik(x-x')+Sik-\t-t') 

-2ifc(x-a;')-84fc3(t-t') -2ife(x-x')-84fc^(t-t') 



And we have the following inequalities on the contours: 



(2.11) 



7i 

72 
73 



x-x' > 0,t-t' < 0, 
x-x' > 0,t-t' > 0, 
X — x' < 0. 



(2.12) 



So, these inequalities imply that the functions {/ij}i are bounded and 
analytic for A; G C such that k belongs to 



/^2 
/^3 



(^2,^2,^3), 
(Di,Di,D4), 
(D3UD4,D3UD4,DiUD2). 



(2.13) 



where {-D„}^ denote four open, pairwisely disjoint subsets of the Rie- 
mann /c— sphere showed in Figure 2. And the sets {-Dnji has the fol- 



J.XU AND E.FAN 




Figure 2. The sets D„, n= 1, . . . , 4, which decompose 
the complex /c— plane. 

lowing properties: 

Di = {k E C|Re/i = ReZ2 > Re/3, Rezi = Re 2:2 > Re^s}, 
D2 = {k E C|Re/i = Re/2 > Re/3, Re^i = Re2;2 < Re^s}, 
Di = {k E C|Re/i = Re/2 < Re/3, Kezi = Iiez2 > Re^s}, 
Di = {k E C|Re/i = Re/2 < Re/3, Rezi = Rez2 < Re2;3}, 

where li{k) and Zi{k) are the diagonal entries of matrices —ikA and 
—Aik^A, respectively. 

In fact, for x = 0, /ii(0,t, A;) has enlarged domain of boundedness: 
{D2 U D4,D2 U D4, Di U Ds), and fi2{0,t,k) has enlarged domain of 
boundedness: {Di U D3,DiU D3, D2U D4). 

2.3. The M„'s. For each n = 1, . . . , 4, define a solution M„(x, t, k) of 
(IA.2p by the following system of integral equations: 



{MnUj{x,t,k) = 6,,+ I [e^-^'^-^''^ *)^W^„(x',t',fc)),„ k e Dn, I, J = 1,2,3. 

(2.14) 
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where Wn is given by f l2.9p with fj, replaced with M„, and the contours 

7j"-, n = 1, ... ,4, i,j = 1,2,3 are defined by 

{7i if Re^j(fc) < Iielj{k) and liezi(k) > 'ReZj(k), 
72 if Reli{k) < Relj{k) and Rezi{k) < Rezj{k), for k G D„. 
73 if Rek{k) > Relj{k) 

(2.15) 

The following proposition ascertains that the M„'s defined in this 
way have the properties required for the formulation of a Riemann- 
Hilbert problem. 

Proposition 2.1. For each n = 1,...,4, the function Mn{x,t,k) is 
well-defined by equation (2J^ for k G l)„ and {x,t) G Q. For any 



fixed point {x,t), Mn is bounded and analytic as a function of k E Dn 
away from a possible discrete set of singularities {kj} at which the 
Fredholm determinant vanishes. Moreover, Mn admits a bounded and 
contious extension to Dn and 

Mn{x,t,k) =I + 0{-), k^oo, keDn. (2.16) 

Proof. The bounedness and analyticity properties are established in 
appendix B in [12] . And substituting the expansion 



M = Mo + — — + -—- + ■■■ , k-^oo. 

into the Lax pair (1A.2I) and comparing the terms of the same order of 
k yield the equation (12.161) . D 

2.4. The jump matrices. We define spectral functions Sn{k), n = 
1, ... ,4, and 

Snik)=MniO,0,k), keDn, 71=1,..., 4. (2.17) 

Let M denote the sectionally analytic function on the Riemann fc— sphere 
which equals M„ for k G Dn. Then M satisfies the jump conditions 

Mn = M„iJm,n, k E Dn H Dm, n,m = l,...,4, n^m, 

(2.18) 
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where the jump matrices Jm,n(yX,t, k) are defined by 

T (—ikx—iik^t)K/Q—lQ \ 

According to the definition of the 7*^, we find that 



(2.19) 



r 



r 



73 73 73 \ 

73 73 73 

72 72 73 / 

/ 73 73 7i \ 

73 73 7i 

V 73 73 73 / 



r 



r 




(2.20) 



2.5. The adjugated eigenfunctions. We will also need the analytic- 
ity and boundedness properties of the minors of the matrices {/ij(x, t, k)}l. 
We recall that the adjugate matrix X"^ of a 3 x 3 matrix X is defined 
by 

/ mn(X) -mi2(X) m.^iX) \ 
X^ = -m2i{X) m22{X) -m2-i{X) , 
V m3i(X) -m;2{X) m33(X) / 

where mij{X) denote the {ij)th minor of X. 

It follows from (JA.21) that the adjugated eigenfunction /i^ satisfies 
the Lax pair 



/^x 



fif-[4ik^A,fi^] = -V^fi^. 



(2.21) 



where V'^ denote the transform of a matrix V. Thus, the eigenfunctions 



{fj'f}l are solutions of the integral equations 



/i • (x, t,k) =1 



Jk{x-x')+4ik'^(t-t')A/ 



{V,'dx + V2')fi 



Ij 



1,2,3. 



(2.22) 



Then we can get the following analyticity and boundedness properties: 






{D,,Ds,D2), 
{D^,D^,D^), 
{DiUD2,DiUD2,D3UD^). 



(2.23) 
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In fact, for x = 0, yU^(0, t, k) has enlarged domain of boundedness: 
{Di U D3, Di U -D3, D2 U D4), and /i^(0, t, A;) has enlarged domain of 
boundedness: {D2 U D4, D2 U D4, Di U Dg). 



2.6. The Jm,n''s computation. Let us define the 3 x 3— matrix value 
spectral functions s{k) and S{k) by 



//3(x, t, k) = ^2(x, t, k)e^-'^''-^"'''^^s{k), 



(2.24a) 



Thus, 



fii{x, t, k) = ^i2{x, t, k)e^-'^''-'''''''^^S{k), (2.24b) 



s{k)= fi3{0,0,k), S{k) = fMi{0,0,k). (2.25) 



And we deduce from the properties of Hj and /xf that s{k) and S'(A;) 
have the following boundedness properties: 

s{k): {DsUD4,DsUD4,DiUD2), 
S{k): {D2UD4,D2UD^,D,UD:i), 
s^{k): {DiUD2,DiUD2,DsUD4), 

S^'ik): (DiUA3,AUD3,^2UD4). 



Moreover, 

Mn{x,t,k) 



1^2 



{x,t,k)e^-'''''"^'''''^^Sn{k), keD„ 



(2.26) 



Proposition 2.2. The Sn can be expressed in terms of the entries of 
s{k) and S{k) as follows: 



Si 



So 



I rn22is) m.2i(s) 
S33 S33 

mi2{s) mii{s) 



S33 





S33 




Sl3 

■S23 



\ 



V S33 / 

/ m22 (^)™-33 (>S')-m32 {s)m23 (5) m2i (s)m33 (5)-m3i (8)m23 (S) 

{sTS^)33 (sTSA)33 

Tni2 (s)m33 (S')-m32 {s)mi3 (5) m-n (s)m33 (S)-m3i {s)mi3 (S) 

WS^3 is''S^)33 

mi2(s)m23(5')-m22(^)mi3(5) mn (s)m23 (S)-m2i {s)mi3 (S) 
\ (3^3^)33 (s^5'*)33 



Sl3 
S23 



S33 / 

(2.27a) 
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S. 



Sll Si2 
521 522 



523 



V ^31 532 (s-5i)33 / 



-S4 



' 5ii Si2 

521 522 

V 531 S32 








"133(5) 






(2.27b) 



Proof. Let 7^° denote the contour (Xo,0) — )■ (x,t) in the (x,t)— plane, 
here Xq > is a constant. We introduce fi3{x,t, k; Xo) as the solution 
of fl2.10p with j = 3 and with the contour 73 replaced by 73 " . Similarly, 
we define Mn{x,t,k;Xo) as the solution of fl2.14p with 73 replaced by 
7^°. We will first derive expression for Sn{k; Xq) = M„(0,0, /c; Xq) in 
terms of S{k) and s{k] Xq) = ^3(0, 0, k; Xo). Then (12:271) will follow by 
taking the limit Xq — )■ 00. 

First, We have the following relations: 

M4x,t,k;Xo) = fii{x,t,k)e^-''^-^^'^''^^Rn{k;Xo), 
Mn{x,t,k;Xo) = /i2(x,t, A;)e(-^^^-4*'='*)'^5„(A;;Xo), (2.28) 

M„(x, t, k; Xo) = fisix, t, k)e^-'^--^^'''^^T^{k; Xo). 

Then we get i?„(fc; Xq) and T„(A;; Xq) are fedined as follows: 

e''''^^M^iO,T,k;Xo), 



Rnik;Xo) 

T^{k;Xo) = e"'^^M^{Xo,0,k;Xo). 
The relations (12.281) imply that 

s{k;Xo) = Sn{k;Xo)T-\k;Xo), 



(2.29a) 
(2.29b) 



Sik) = Snik;Xo)R-\k;Xo). 

(2.30) 
These equations constitute a matrix factorization problem which, given 
{s, S} can be solved for the {Rn, Sn, Tn}. Indeed, the integral equations 
(I2.14P together with the definitions of {Rn, Sn, Tn} imply that 



(2.31) 



iRnik;Xo))^,=0 if 75 = 71, 
{Sn{k-Xo))i,=Q if 7." =72, 
(T„(fc;Xo))., = if 75 = 73. 

It follows that (I2.30p are 18 scalar equations for 18 unknowns. By 
computing the explicit solution of this algebraic system, we find that 
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{Sn{k; Xo)}f are given by the equation obtained from f l2.27p by replac- 
ing {Sn{k),s{k)} with {Sn{k; Xo),s{k; Xq)}. taking Xq — )■ cxd in this 
equation, we arrive at (12.271) . D 

2.7. The global relation. The spectral functions S{k) and s{k) are 
not independent but satisfy an important relation. Indeed, it follows 
from fl^:^ that 

/ii(x, t, k)e^~'''''-^'''''^^S-\k)s{k) = /i3(x, t, k), k G {D3UD4, D3UD4, D1UD2 

(2.32) 
Since /ii(0,T, A;) = I, evaluation at (0,T) yields the following global 
relation: 

S-\k)s{k) = e^''''^^c{T,k), k e {D3UD^,D3UD^,DiUD2). (2.33) 

where c{T, k) = /i3(0, T, k). 

2.8. The residue conditions. Since /i2 is an entire function, it follows 
from (I2.26P that M can only have sigularities at the points where the 
S'^s have singularities. We infer from the explicit formulas 02.27P that 
the possible singularities of M are as follows: 

• [M]i could have poles in Di U D2 at the zeros of s^^^ik); 

• [M]i could have poles in D2 at the zeros of {s^ S^)^^{k)] 

• [M]2 could have poles in Di U D2 at the zeros of S33(A;); 

• [M]2 could have poles in D2 at the zeros of (s'^S''^)33(A;); 

• [M]3 could have poles in D^ at the zeros of (5'"^s'^)33(/c); 

• [M]3 could have poles in D^ U D^ at the zeros of m33(s)(A;); 

We denote the above possible zeros by {kj}i and assume they satisfy 
the following assumption. 

Assumption 2.3. We assume that 

• s^^{k) has riQ possible simple zeros in Di denoted by {%}""; 

• ■S33(fc) has ni—UQ possible simple zeros in D2 denoted by {kj}^l^i; 

• (s^5'^)33(/c) has 712 — ni possible simple zeros in D2 denoted by 

l^lni+l/ 
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• (<5'^s'^)33(A;) has ns — 77-2 possible simple zeros in D3 denoted by 

• 77233(5) (A;) has 774 — 773 possible simple zeros in D3 denoted by 

i^iJna+l'' 

• 77733(5) (A;) has 775 — 774 possible simple zeros in D^ denoted by 

l%ln4+l'' 

• "^33 (s) (A;) has N — n^ possible simple zeros in D4 denoted by 

i^lns + l/ 

and that none of these zeros coincide. Moreover, we assume that none 
of these functions have zeros on the boundaries of the Dn 's. 

We determine the residue conditions at these zeros in the following: 



Res,=k,[Mh = .'^^l^^^^)ly ''^''\Mik,)],, 1 < J < 770, k, E D, 

S33['^j)Sl3[t^j) 

Res - \M] = "^12(s)(fcj)m33(S)(fcj)-m32(g)(fcj-)mi3(5)(fcj) g^^(fc.)r^^/^^ N-i 



Proposition 2.4. Let {M„}^ be the eigenf unctions defined by 1^2. 14\ ) 
and assume that the set {kj}f of singularitues are as the above as- 
sumption. Then the following residue conditions hold: 

Res,=,^[Ml = T^^^^^^^e^-(..)[M(A:,.)]3, 1 < J < ^o, k, G D, 

S33{kj)S23{kj) 

(2.34a) 

(2.34b) 

kj)]3 

77i + 1 < j < 772, A;j G -D2, 

(2.34c) 

ReSk-k [Ml = "^^i(^)fe)"^33(.S')(fcj)-m.3i(g)(fcj)m23(g)(fcj) ^e32(fcj)r;^^/^g \1 
~j'- {s'^ S^)33{kj)sx3{kj) ^'^ 

rii + l < j < n2,kj e D2, 

(2.34d) 

(2.34e) 
Resk=k[M]3 = ^^-ii2M^^e^«('=^)[M(A;.)]i - -^wr^^VwTe'''^''-'^[^(^7)]2 

K Kji JO m33(s){kj)m23(s)(kj) L V J/JJ- m33(s){kj)m23(s)(kj) L V JJl^ 

(2.34f) 
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where f = -J^, o,nd 6ij is defined by 

9ij{x,t,k) = {k- lj)x + {zi- Zj)t, i,j = 1,2,3. (2.35) 

that implies that 

% = 0, z,i = 1, 2; ^13 = ^23 = -O32 = -031 = -2ikx - 8tkH. 

Proof. We will prove (EJHD, dOic]) . (l231e|) . dmij), the other condi- 
tions follow by similar arguments. Equation ( I2.26P implies the relation 

Ml = //ae^-*'"-'''^'*)^^!, (2.36a) 

M2 = 1^26^-'''^-^"'''^^ S2. (2.36b) 

M3 = fX2e^-''^-^"'''^^S3, (2.36c) 

M4 = fi2e^-'''''-^"'''^^Si, (2.36d) 

In view of the expressions for 5*1 and S2 given in (I2.27p . the three 
columns of 02.36ap read: 

m, = Hi^^^^ + N^e^-^^^i^, (2.37a) 

S33 S33 

[M1J2 = [^2jie h [/U2J2 , (2.37b) 

■533 ■S33 

[Mi]3 = [/i2]ie''"si3 + [fi2]2e'''s23 + [fi2]3S33- (2.37c) 

while the three columns of (12.36bp read: 

FA/f 1 — \,, 1 "^22 (s)m33 (■5')-m32 (s)m23 (S) 
[M2J1 - [^2jl (sT5A)33 

+ [/X2] 2 "^^ (s)m33 (S)-rn32 (s)mi3 (5) ^^21 (2.38a) 

iF,, 1 mi2(s)m.23(5')-m.22(^)mi3(S) ^g,^-, 
+ [H3 (s'J'5^)33 ^ 

[M2]2 = [/X2] 1 ""'' ^'^""'\[VsT)ll ^'^"'^ ^^^ e^" 

+ [//2] 2 ""^'^"^'(^^51^33 ^'^"" ^^^ (2.38b) 

+ [H3 (sJ'5^)33 ^ 

[M2]3 = [/i2]isi3e'^^ + [fi2hs23e''' + [/i2]3S33. (2.38c) 

and the three columns of (I2.36cp read: 

[M3]i = [Hisii + [l^2]2S2ie''' + [/x2]3S3ie'^^ , (2.39a) 
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-,^12 I [,, 1 „ I [,, 1 o „«32 



[Ms]2 = [^l2]lSl2e''' + [fl2]2S22 + N3532e'^='% (2.39b) 

mh = [f,2]i , J'' e'-^ + [iX2]2 ,J'' e'^' + [f^2h ,J'' . (2.39c) 



(5^5^)33 '"^^'^S^s^U ""^^'^S^s 
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the three columns of fl2.36dp read: 

[Mji = [fi2]isu + [fi2]2S2ie''' + [t^2]3S3ie''\ (2.40a) 

[M4]2 = NlSl2e'" + [H2S22 + [^35326^^% (2.40b) 

[M.h = [/^2]3-^. (2.40c) 

We first suppose that kj G -Di is a simple zero of s^si^k). Solving (I2.37cp 
for [/i2]2 and substituting the result in to (I2.37al) . we find 

[Miji = e ^'[Mijs H [^2j2 e ^'[^2j3- 

■S33'S23 ■§23 ■§23 

Taking the residue of this equation at kj, we find the condition (I2.34ap 
in the case when kj G Di. Similarly, Solving (I2.38cp for [^^2)2 and 
substituting the result in to (I2.38ap . we find 

(5 "J J33'S23 -523 S23 

Taking the residue of this equation at kj, we find the condition (I2.34cp 
in the case when kj G D2. 

In order to prove (I2.34ej) . we solve 02.39aj) and (]2.39bp for [/i2]i and 
[/i2]3, then substituting the result into (I2.39cp . we find 

[M3]3 = f"X ^"f. e^-[M3],+ . J^^;." - ^"^" e^^^lMM^M- 

(S^S^)33m23(s) {S^S^)33{kj)m23{s) ^23(5) 

Taking the residue of this equation at kj, we find the condition (I2.34ep 
in the case when kj G D3. Similarly, solving (I2.40ap and (]2.40bp for 
[/i2]i and [/i2]3, then substituting the result into O2.40cp . we find 

mh = ''' , , e'-^[M,], '" ,/ ''[M,]2 V^e^-H2. 

"^33(s)"^23(s) ^33(5)77123(5) m23(s) 

Taking the residue of this equation at kj, we find the condition fl2.34fp 
in the case when kj & D^. D 
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3. The Riemann-Hilbert problem 

The sectionally analytic function M{x, t, k) defined in section 2 sat- 
isfies a Riemann-Hilbert problem which can be formulated in terms of 
the initial and boundary values of u{x,t). By solving this Riemann- 
Hilbert problem, the solution of fll.Sp fthen (11.31) ) can be recovered for 
all values of x, t. 

Theorem 3.1. Suppose that u{x,t) is a solution of U.5\) in the half- 
line domain Q with sufficient smoothness and decays as x -^ oo. Then 
u{x, t) can be reconstructed from the initial value {uo{x)} and boundary 
values {gQ^t) , gi(t) , g2{t)} defined as follows, 

Uo{x) =u{x,0), go{t) = u{0,t), gi{t) = u^{0,t), g2{t) = u^x{0,t). 

(3.1) 

Use the initial and boundary data to define the jump matrices Jm.n{x, t, k) 
as well as the spectral s{k) and S{k) by equation 1^2.24^ . Assume that 
the possible zeros {%}^ of the functions s^^{k)^ {s^ S'^)^3{k)^ {S^ s'^)^3{k) 
and m3^{s){k) are as in assumption 2.3. 

Then the solution {u{x,t)} is given by 

u{x,t) =2i\im{kM{x,t,k))i3. (3.2) 

where M{x, t, k) satisfies the following 3x3 matrix Riemann-Hilbert 
problem: 

• M is sectionally meromorphic on the Riemann k~sphere with 
jumps across the contours Dn^Dm, n, m = 1, ■ ■ ■ ,4, see Figure 
2. 

• Across the contours Dn fl Dm, M satisfies the jump condition 

Mn{x, t, k) = Mm{x, t, k)Jm,n{x, t,k), k E Dn H Dm, 71,171 = 1, 2, 3, 4. 

(3.3) 

• M{x,t,k) =1 + 0(1), k^oo. 

• The residue condition of M is showed in Proposition\2.4- 
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Proof. It only remains to prove (13.21) and this equation follows from 
the large k asymptotics of the eigenfunctions, see the appendix \M D 



4. NON-LINEARIZABLE BOUNDARY CONDITIONS 

A major difficulty of initial-boundary value problems is that some of 
the boundary values are unkown for a well-posed problem. All bound- 
ary values are needed for the definition of S{k), and hence for the for- 
mulation of the Riemann-Hilbert problem. Our main result expresses 
the spectral function S{k) in terms of the prescribed boundary data 
and the initial data via the solution of a system of nonlinear integral 
equations. 



4.1. Asymptotics. An analysis of flA.2p shows that the eigenfunctions 
{fj'j}l have the following asymptotics as A; — )■ oo (see the appendix R|) : 



LljiXj Z J th j 



i + i 



+^ 



+ fc3 



where 



/ 



V 



-(x,t) 
2 J{xj,tj) 






A 

A2 



J(x..U) ^1 



-{x,t) 

l{Xj,tj) 

i r{a;,t) ^ 



' {Xj,tj) 



h^ 



^4£l)(^+-i) -ij; 



_1 f(a:,t) 
4 J( 






{x,t) 

(Xj,tj) 



Vi 



( u^') 

' Pll 

V 4? 



P12 

(3) 
/^22 

(3) 
^^32 



(3) 

(3) 
/^23 

1^1 I 



-j/,':l,(i+''.) 



(2) 
^32 



(Xj,tj) 



+ o{h 



21 
2i 



\ 



J(xj,tj\ 



A 



(2) 
/^13 

(2) 
/^23 

^{x,t) 



\ 



r(x,t) 
Jix„t,)V3 ) 



(4.1a) 



A = 

Ai 



Ittpdx + {uUxx + UxxU — u^Ux + 6\u\'^)dt 
-u^dx + {uuxx + UxxU — {uxY + Q\u\'^v?')dt 



A2 = -U^dx + [uUxx + MxxM 



'U.T 



+ 6|ti|^w^)rft 



(4.2a) 
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{x,t) 



2"' J(x,.U) ^ 



4^a; 



(2) 

(2) 
1^23 



2"' J(xj,t/) ^ ^ 4"^ 



u.,.. 



(4.2b) 



ui = Ai /(^^^' j^) A2 + {uu^)dx + {uut - 2\u\'^{uu^ - u^u) - {u^^Ux - UxUxx))dt 
= r(^.*) 'X r^^''*^ Ai + u^ 

_ r(a;,t) A r(2^-,t) a , -2 
^2 - J(x,,t^)^J(r.^.^^2+M 

^{x,t) 



{x,t) /^ ^(a-,t) 
i '*i ) 
1^2 = A2 /(If ''(,) Ai + {uUx)dx + {uUt - 2\u\'^{uU^ - U^u) - {UxxUx - UxUxx))dt, 
{x,t) ^x2 

(4.2c) 

(3) ,.(3) ,.(3) _^ „(3) 



V3 



d[-Ul!xl^ 



ihpl 



and in the following we just use /uia , /i23 , /^-ji and fi^ , so we only 
compute these functions 






(3) 
/^31 

(3) 
/^32 



(tt/i 



(2)1 



{2)^ 
21 y 

(2). 



1, 



2.,.^ ^"+^/i^7)-i(^l,/X«" 



«xMV) + i|MpU + |tia;x 



^(m^(2)i2 +^x//^-' 



i(^./^(^) 



12 +U^^^22) + j\u\'^U + ^U^ 



(4.2d) 



From the global relation (I2.33p and replacing T by t, we find 



^2(0,t,A;)e-^*'='*^s(A;) = c{t,k), ke {D3U D4, D3U D4, DiU D2). 

(4.3) 
We define functions {$i3(t, k), $23(t, k), $33(t, A;)} and {cj(t, A;)}^ by: 



fi2{0,t,k) 



/ <^nit,k) $i2(t,A;) $i3(t,A;) \ 

$21 (t, A:) <l>22(t,A;) $23(t,A;) 
V $31 (t, A:) $32(t,A;) $33(t,A;) / 



[cit, 



S33{k) 



ci{t,k) \ 

C2{t,k) 
Cs{t,k) J 

(4.4) 



we can write the (13) and (23) entries of the global relation as 



$n(t, fc)e-«^'='*^+<l'i2(t, k)e-''''''^+^,,{t, k) = c,{t, k), ke D^UD^, 
S33 S33 

(4.5a) 
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$21 (t, k)e-^'^"'^+^22{t, k)e-^'^''^+^23it, k) = C2(t, A;), ke DiUZ^s, 



■S33 



■533 



(4.5b) 

The functions {cj{t,k)}\ are analytic and bounded in Di U D2 away 
from the possible zeros of S33(A;) and of order O(^) as A; — t- 00. 
From the asymptotic of /ij(x, t, A;) in (14. lap we have 



and 



S23ik) 
\ S33{k) J 



$,3(t,A:) 



/0\ 





u 



(0,0) \ 



+ 



2ik 



m(0,0) 
rim 

/(oo,0) 



2 r^°'°^ A / 



+o(t,; 



(4.6) 



k k"^ k^ 



«'f)- 



(4.7a) 



$33(t,A;) = l+^|i^ + ^|fi+0(-), k ^ 00, k&D.VJD^. (4.7b) 
where 



* 3^(^) = iMtY. 



^>) = \9.{tY-\9lJtS^ 



^'^{t) = yi^f^ + \gl^^ + l\u\-gl + lgl 



(0,0) 



.(1) 



^"^ify 



■ ri^,t) ^ 



.(2) 



^!^w = i;r 



{x,t) 

(Xj ,tj) 



m- 



Here the definition of $j3(t, A;) can be found in the appendix IM 

In particular, we find the following expressions for the boudary val- 
ues: 

g^ = 2z<^^^it), (4.8a) 

^f = 2^^o^$g(t)+4$f3)(t), (4.8b) 

g^ = -2\go\'g^ + 2^^f $g(t) + Ag^^f^it) - Sz^f^it). (4.8c) 

We will also need the asymptotic of Cj{t, k), 

Lemma 4.1. The global relation ( [/^.5] ) implies that the large k behavior 
of Cj{t, k) satisfies 



Cj{t,k) 



<^S(o , <^s?w , *;?w 



+ 



k k 

Proof. See the appendix [ 



fO(-), A;^oo,A;eDi. (4.9) 



D 
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4.2. The Dirichlet and Neumann problems. We can now derive 
effective cliaracterizations of spectral function S{k) for the Dirichlet {go 
prescribed), the first Neumann {gi prescribed), and the second Neu- 
mann [g2 prescribed) problems. 

Define ahj a = e^~ and let {n.j{t, k), Ilj{t, k), Ilj{t, k)}^ denote the 
following combinations formed from {<l>j3(t, A;)}'^: 

Uj{t, k) = ^j3it, k) + a$j3(t, ak) + a'^^j3{t, a^k), j = 1, 2, 3, 
Ilj{t, k) = ^j3{t, k) + a^^j3{t, ak) + a$j-3(t, c^k), j = 1, 2, 3, 
n,(t, k) = <l>,3(t, A:) + $,-3(t, ak) + <l>,-3(t, a^k), j = 1, 2, 3. 

(4.10) 

Andleti?(A;) = $n|f + $i2|f. 

Let Di = D[\JD'[ where D\ = Din{Rek > 0} and D'^ = Din{Rek < 
0}. Similarly, let D^ = D'^U DI where D'^ = D^r] {Rek > 0} and 
Dl = DiH {Rek < 0}. 

Theorem 4.2. LetT < oo. Letuo{x),u > 0, be a function of Schwartz 
class. 

For the Dirichlet problem it is assumed that the function gQ(t),0 < 
t < T, has sufficient smoothness and is compatible with Uo{x) at x = 
t = 0. 

For the first Neumann problem it is assumed that the function gi (t) , < 
t < T , has sufficient smoothness and is compatible with Uq{x) at x = 
t = 0. 

Similarly, for the second Neumann problem it is assumed that the 
function g2{t),0 <t<T, has sufficient smoothness and is compatible 
with Uq{x) at X = t = 0. 

Suppose that s^slk) has a finite number of simple zeros in Di. 

Then the spectral function S{k) is given by 



Sik) 



/ A{k) B{k) e^^'^^'^Cik) \ 

D{k) E{k) e^'^'^F{k) 

\ e-^"''^G{k) e-^'^'^H{k) I{k) j 



(4.11) 
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where 

A{k) = ^22{k)^zm - <I>23(A:)$32(A;) B{k) = $i3(A:)$22(A:) - $12 (A:) $33 (A;) 

C{k) = $i2(A;)$23(A;) - ^n{m22{k) D{k) = $23 (A;) $31 (A:) - $21 (A;) $33 (A:) 

^(A:) = $n(A;)$33(A:) - $i3(A;)$3i(A:) F{k) = <l>2i(A:)<l>i3(A;) - $n(A:)<l>23(A;) 

^(A;) = <l>2i(A;)<l>32(A;) - $22 (A;) $31 (A;) H{k) = $i2(A;)<l'3i(A;) - <l>n(A:)<l'32(A;) 
/(A;) = <l>ii(A;)$22(A:)-$i2(A;)$2i(A:) 

and the complex-value functions {$i3(t, A;)}f^^ satisfy the following sys- 
tem of integral equations: 

$i3(t, A;) = /,* e-8^^'(*-*') [i2tk\go\' + (Mi " gi9o))^i3 

+^2$23 + {Ak^go + 2ikgi - 4|^o|'^o - ^2)^33] {t\ k)dt' 

(4.12a) 

$23(t, A;) = Jl ^^^'^^~^) [{2ik\g,\' - {g.g, - ^i^o))$i3 

+91^23 + (4A:2^o + '^ikgi - 4|(7oP^o - ^2)^33] it\ k)dt' 

(4.12b) 

<l>33(t, A:) = 1 + /q [(-4A;2^o + ^^kg, + 4\go\'' + ^2)<5i3 

+ (-4A;2^o + ^ikgi + 4|^oP + ^2)^23 + -4iA;|^oP<^33] {f, k)dt' 

(4.12c) 

and {(^iiit, A;)}f^^, {$^2(^1 ^)}f=i satisfy the following system of integral 

equations: 

$ii(t, A;) = 1 + jl [(2iA;|t/oP + (Mi - ^i^o))<^ii 

+gl^2i + (4A;2^o + 2ikgi - 4|^o|'^o - ^2)^31] it', k)dt' 

(4.13a) 

<^2i(t, A;) = /p [(2^A;|^oP - (Mi - 9i9o))^ii 

+gl^2i + (4A:2^o + 22A;^i - 4|(/o|'^o - ^2)^31] it\ k)dt' 

(4.13b) 

$33(t, A:) = /J eS*'^'^*-*') [(-4A:2^o + 2tkg^ + A\go\^ + ^2)^11 

+ (-4A;2^o + 2«A;^i + 4|^oP + ^2)^21 + -4iA;|^oP'^3i] (f, k)dt' 

(4.13c) 

$i2(t, A;) = jl [(2«A;|^oP + (Mi - 9i9q))^i2 

+gl^22 + (4A;2^o + 2ikgi - 4|^o|'^o - ^2)^32] (f, A;)dt' 

(4.14a) 
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$22(t, A;) = 1 + /o* [{2ik\go\'^ - {goQi - gigo))<^i2 

+gi<l>22 + (4A;2^o + 2ikgi - 4|^o|'^o - ^2)^32] (f, k)dt' 

(4.14b) 

$32(t, A:) = /J e8*'^'(*-*') [(-4P^o + 2tkg^ + 4|(7oP + ^2)^12 

+ {-Ak^go + 2ikgi + 4|^oP + ^2)^22 + -4iA;|^oP<^32] (f, A;)rft' 

(4.14c) 

(i) For the Dirichlet problem, the unknown Neumann boundary 
values gi{t) and g2{t) are given by 



9iit) 



^ /,^^ n3(t, k)dk + A /,^^ \ku,{t, k)-'-^ 



2i 



dk 



,-8ik-h\(„,2 



[(a^ - a)R{ak) + (a - a^)i?(a2A;)](iA; 



+4 {(1 - «^) E.^ez.; +(1 - «) T.u,eD'i } k,e~'^'^'Res,^Rik). 

(4.15a) 



an(i 



92{t)= goitf-tlan.. k'Mt^k)- 



+;; JaDs ^ 



2„-8ifc^t 



3fcgo(t) 
2j 



t/A; 



[(1 - a)R{ak) + (1 - a2)i?(a2^)] ^^A; 



TT JdD^ Lv / V / ' V / V /J 

-8z {(1 - «) E.^eDi +(1 - «') E,,eD';}k]e~''''^'Res,^R{k) 



4go(t) 



/,^^^ m3(t, A;)rffc + ^ /,^^ n3(t, A;)rfL 



(4.15b) 

(ii) For the first Neumann problem, the unknown boundary values 
go{t) and (72 (t) are given by 



g^(t) = i Jgj^^ ni(t, k)dk - i /g^^ e-s*'^'* [(a - alpha^)R{ak) + (a^ - a)i?(a2A;)] dA; 

R{k), 
(4.16a) 



-2^ {(1 - «) E,,en[ +(1 - «') E.,eD'/ } e-''''''Res,^R{k), 



and 



92{t) 



9l^t) - -. Ian, {k'Mt, k) - ^^ /,^^ l^t, l)dl) dk 

+ ! !dD, A;2e-8^'^'* [(1 - a')R{ak) + (1 - a)R{a^k)\ dk 

-8z {(1 - a) Y.u,eD[ +(1 " «') E.,,z.;' } ky'^'^' Res,^R{k) 



/,. m3(t, A;)rffc + ^ / n3(t, A;)rfL 



TTJ JdDi 



(4.16b) 
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(iii) For the second Neumann problem, the unknown boundary val- 
ues go{t) and gi{t) are given by 

g^{t) = i /g^^ Ili(t, k)dk - i /g^^ e-s*'^'* [{a - alpha^)R{ak) + {a^ - a)R{a^k)] dk 
+2^ ((1 - «) Y.u,eD[ +(l'- «') E,,ez.;' } e-'^'^' Res,^R{k) , 



(4.17a) 



and 

_ -290 jt) 



9i{t) = ^ Jan., Mt, k)dk + ^ /,^^ kn,{t, k)dk 

-Vi IdD3 ke-^'^'' [{a^ - l)R{ak) + (a - l)R{a^k)] dk 

+4{(i-«)E.,ei.;+(i-«')E. 



k,eD';fk,e-'^'^'Res,^R{k). 
(4.17b) 

Proof. The representations (14. lip follow from the relation S{k) = e^**^ "^//^(O, T, A;)^. 
And the system (I4.12p is the direct result of the Volteral integral equa- 
tions of /X2(0, t, k). 

(i) In order to derive (14. 15ap we note that equation f l4.8bp ex- 
presses gi in terms of $33 and "^13 . Furthermore, equation 
(14. 7p and Cauchy theorem imply 



2m 



^(l) 



^'33' it) 



2m 



3 

and 

Thus, 



2 / [$33(t, k) - l]dk = / [<l>33(t, A;) - l\dk 
'dD2 JdDi 



dD2 L 



k^n{t) 



9o{t) 
2i 



dk 



dDi 



k^n{t) 



9o{t) 
2i 



dk. 



Z7r$g(i) = - {I9D2 + Ld,) ['^33(t, k) - l]dk = (/,^^ + /,^J [$33(t, k) - l]dk 

= /a^3[$33(t, k) - l]dk + a/g^J$33(t, k) - \\dk + a" /ao3[$33(t, A:) - V^dk 
= IdD3'^^(^'k)dk. 

(4.18) 

Similarly, 



^^^rJit) = ilen, + IoDj k^^.it) 



go(t) 



2i 
S0(t) 



dk 



l9D.+^'IaD[+^IaD';) k<^rs{t)-'-^ dk + m (4.19) 



Ld. kMt,k) 



3go(f) 

2i 



dk + I{t). 
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where I{t) is defined by 



m 



a 



a] 



dD\ 



dD'l 



k^nit) 



2i 



dk 



The last step involves using the global relation to compute I{t) 



m 



a 



h 



do: 



'I -a) 



IdD';) kci{t,k) 



9o{t) 
2i 



dk 



((1 - a') /,^, +(1 - a) /,^„) ke-'^>^'^R{k)dk 



(4.20) 



Using the asymptotic (14. 9 p and Cauchy theorem to compute 
the first term on the right-hand side of equation (I4.20p and 
using the transformation k —^ ak and k — )■ a'^k in the second 
term on the right-hand side of (14.201) . we find 



I(t) = -in<!?\y{t) - Jg^^ ke-^'''' [{a^ - a)R{ak) + (a - a^)R{a'^h)] dk 

+2^^ {(1 - «') E.,eD; +(1 - «) Y.,,^o'l]ke-'''^'Res,^R{k). 

(4.21) 
Equations fj4T9|) and fj42TD imply 



"^13^ = 2^/0^3 kn,{t,k)- 



IdDs 

u3 



3go{t) 

2i 



dk 



-h ho, ke-^'^' [(a^ - a)R{ak) + (a - a^)R{a'^k)] dk 
{(1 - «'') E.,eD; +(1 - «) i:u,eD'l}k,^-'''^'Resu,R{k)- 



This equation together with ( ]4.8b[) and (14.18P yields ( 14.15ap . 

In order to derive fl4.15bp . we note that (I4.8cp expresses g2 
in terms of $13 , $33 and $33 . Equation fl4.15bp follows from 
the expression (I4.18p for $33 and the following formulas: 



.(2) 



^tiit) 



m 



kll-idk, 



(4.22a) 



dDs 



M, 



3fcgo(f) 
2i 



dk 



*i3W = 2^Urni(t,fc) 

-2^ IdD, k'^e-^'^'' [(1 - a)R{ak) + (1 - a'^)R{a'^k)] dk 
1 - «) E.,eD; +(1 - «') j:,,eD';}k]e~'^''^'Res,^Rik). 



(4.22b) 
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(ii) In order to derive the representations f l4.16p relevant for the 
first Neumann problem, we use (14 .Sp together with fl4.18p . 
fl4.22al) and the following formulas: 

-"ti IdDs e"^^'''* [(a - a^)R{ak) + {a^ - a)R{a^k)] dk (4.23a) 

$g(t) = — / ktiidk, (4.23b) 

7r« JdD-, 



'dDs 

(2) 



dk 



^rJit) = ^J,j,^[k'U,{t,k)-3^\^, 

-h IdDs k^e-''"''' [(1 - a^)R{ak) + (1 - a)R{a^k)] dk (4.23c) 

{(1 - «'') E,,eD; +(1 - «) j:,,eD'^}k]e-'''^'Res,^R{k). 

(iii) In order to derive the representations (I4.17P relevant for the 
second Neumann problem, we use (14. 8 p together with (I4.18P 
and the following formulas: 

-h IdDs e"''"'* [(« - l)R{ak) + (a2 _ l)R{a''k)] dk (4.24a) 

{(1 - «') E,,eD[ +(1 - «) E,,eD';}e~'"'^'ReSkMk)- 

^?ht) = ^Jan,km^k)dk 

-h IdDs ^e-''*^'* [(«' - 1)^(«^) + (a - l)R{a^k)] dk (4.24b) 



D 



4.3. Effective characterizations. Substituting into the system (14.121) 
the expressions 



$.,= 


go = egm + £^5-02 H , 

gi = egii + e^gi2 H , 

g2 = eg2i + £^^22 + • • • , 


= 1 


2 


3. 


(4.25a) 

(4.25b) 
(4.25c) 
(4.25d) 
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where e > is a small parameter, we find that the terms of 0(1) give 
<|)(J = $^J = and $^3^ = 1. Moreover, the terms of 0(e) give $^3^ = 
and 

0(e) : <l>^Sit,k) = [\-''''^'~''\4k'goi + 2tkgn- g2iW,k)dt', 
Jo 

(4.26) 

From the above equation (14.261) we can get 



Ilf\t,k) = 12k- I e 
'0 



-Sik'^it-t') 



9oi{t')dt', 



(4.27a) 



Il^^\t,k) = 6ik [ e-^*'='(*-*')(?ii(t')rft', (4.27b) 

Jo 



Ilf\t,k) = -3 / e 
'0 



-'"'''^'-''^gii{t')dt', 



(4.27c) 



The Dirichlet problem can now be solved perturbatively as follows: 
assuming for simplicity that S3s{k) has no zeros and expanding (I4.15ap 
and (I4.15bp . we find 



^11= ^U kn^^\t,k)-'-MMl 



dk 



Vi IdDs ke-^'^'^lia"^ - a)si3i{ak) + {a - a'^)sui{a'^k)]dk 

(4.28a) 



5'21 



.2tt(1), 



-i/,,^ eu[^>it,k) 



3fcgoi ft) 
2i 



dk 



+ t IdD, k'^-'"'"' [(1 - a)si3i(aA;) + (1 - a^)sni{a^k)] dk 



(4.28b) 

Using equation (I4.27ap to determine 11-^ , we can determine gu,g2i 
from (I4.28p . then $^3 can be found from (I4.26p . And these arguments 
can be extended to higher orders and also can be extended to the sys- 
tems (I4.13ap and (I4.14ap . thus yields a constructive scheme for com- 
puting S{k) to all orders. 

Similarly, these arguments also can be used to the first Neumann 
problem and the second Neumann problem. That is to say, in all cases, 
the system can be solved perturbatively to all orders. 
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Appendix A. The asymptotic behavior of the functions 

{Hj{x,t,k)}l 

We denote some symbols as follows: 

A = I '-^ °_^ I , (A.la) 



' -U 



V,^'^ = 2z 



-U 



2 



' ' -f/ ; 1^ _[/^.^. ; ^ ^ ^ ^ 1^ 

(A.lb) 

where I2x2 = I 1 and U = (u,u). 

' 1 ' 



From the Lax pair of /i 



//^ + [ikA, fj] = Vifi, 
lit + [4zA;3A, /^] = \/2/x. 



(A.2) 



Suppose that 



Mx,t,A:) = Do + ^ + ^ + ^ + --- • (A.3) 



We substitute the equation (1A.3P into the Lax pair (lA.2p , and compare 
the order of A;, we find that: 

0{k): [2A,Do] = 0, 

0(1): Do. + [^A,Di]=yiDo, ,^^. 



0(P) : 


[At A, Do] = 0, 


0{k') : 


;4iA,Di] =y2^^)L)o, 


0{k^) : 


4tA,D2]=vi^^D, + vi'^Do, 


0(1): 


Dot + [42A, D3] = ^2 '^^2 + ^2'^Di + VJ^'^Do 


0{k-') : 


Du + [4iA, D4] = V^'^'^D^ + y2^^^D4 + vi'^'^Di 


0(A;-2) : 


D2t + [42A, D5] = vi'^^D^ + ^2(^^1)3 + ^^°^^2 



(A.4b) 



And we denote the Di by D; 
Then, from 0(A;'^),we have 




D 
D 



(0 

(0 
33 



J = 1,2. 



D 



(0) 

i3 



0, 



^3i 



(A.5) 



0{k'^), we get 



4i 



i(i) 



2D 



(1) 
i3 



2D^,/ 



-?7D 



(0) 
2x2 



t/^D(°) 



'33 



(A.6a) 



this imphes that 



D 



(1) 
i3 



^t/^D 



(0) 

33 



3j ~ 2^ 2x2- 



(A.6b) 



0(A;), we find 



4i 



1(2) 



2D 



'3i 

/7^D« 

-t^^S2 



this imphes that 



(2) 

i3 



2Di:.^ 



i3 



-f/D 



+ 2i 



-u d(°) 






(A.7a) 



^j3 
^3j 



2'-^-'^2x2 4L^xi^2x2- 



(A.7b) 
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0(1), we have 







-4|m| 



f/^Dl°) \ / ULDf^ 



XX 33 



Ud'^l / I -U^.D^^l 



2x2 " / \ '-^3::x-L^2x2 



^3-^^2x2 



(A.8a) 



this iinphes that 



3j ~ 2 2x2 4*^ i: -^2x2 T^ 4 I "I '-^ -^2x2 ^^ 8 ^^2x2- 



(A.8b) 



0(A; ^), we get 



-f/DS2 -f/I^f ; A -^^^?x2 - 2\u\^D^S UxD^ - 2|-P^:^'^ 



-4|w| 



+ {UUx — M^iW 



this imphes that 





_r/ n« _f7 n« 



(A.9a) 



^2x2t = fl^^^xx + t/xx^ - U^U^x + 6|lz|2[/^f/}DS2 

^33t = -i{uu^x + U^xU - tx^M^ + Q\u\^}D^^^. 

n(3) _ _±/7 n{0) _ifTr)^3) _ 1 rr r,(2) , iL,|2r7r)(l) , i/y r)(l) _ iL.pTy r)(0) 
3i ~ I6^t-'^2x2 2 2x2 4'-^a:-t^2x2 "•" 4 1 "I "-^ -'^2x2 "•" 8 ^^ 2x2 gl I ^a;-'^2x2- 

(A.9b) 
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0{k^'^), we get 



(A. 10a) 
this implies that 

Asst = -i{uuxx + UxxU - UxUx + Q\u\'^}D)^^ - \{\u\'^)tDf^ . 

(A.lOb) 

Also, from the x— part of the Lax pair, we have the following equations 

DfL = ^. Dfl = 0. (A.Ua) 



-^2x2x ~ 2 ^-^^2x2 
-'^2x2:r ~ 2^ '-^-'^2x2 i^ '-^a:-t^2x2 

D?l = ^\u\'D^^-\{\u\%Df 



(A.llb) 
(A.llc) 



4VI"'l ;a;-t^33 • 

Then from the integral contours 7j, we can get 

^^2 = 12x2, Df;^ = l. (A.12) 

Appendix B. The asymptotic behavior of Cj{t,k) 
Let 

fi2{0,t,k) = 

The global relation shows that 

$2x2^e-«^^^* + <|.,3 = c,. (B.i: 

■S33 
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And from equation 



Ht+ [4iPA,/i] = V2H. 



we get 



$2x2 $i3 



$ 



3j 



$ 



+ 4zP 



33 



2$ 

-2<l>3, 



JS 



4A;2 



+2ik 



?7^f/$2x2 + t/J<^3, f/^^$,3 + f/J$33 



^x'^2x2 — -^1"! ^3j 



2|m|2$ 

— ^xx'^2x2 —Uxx^j3 



-[/,<I>,3-2|m|2$33 



f/^$3, 
-f/$2x2 

-4bP 



U^<!> 



33 



-f/$ 



i3 



f/^$3, 



2x2 



+ (MM^r - UxU) 



Cr3$2x2 Cr3$j3 




fB.2) 



Prom the second column of the equation (1B.2P we get 

$,3t + 8te^,3 = 4A;2[/^<|.33 + 2lk{U^U^,s + [/J$33) 

-4|M|2f/^<|>33 - ?7j^.$33 + (WM^ - UxU)a3^j3 
$33* = -4P[/$,3 + 2lk{U,<l>j3 - 2|u|2$33) + 4|M|2f/$,3 + [/,,$,.3. 

(B.3) 
Suppose 



U^^ 



33 



-f/$ 



i3 



$33 / k k"^ k k"^ 

(B.4) 

where the coefficients ai(t) and A(0? ^ > 0, are independent of k. 
To determine these coefficients,we substitute the above equation into 
equation ( IB. 31) and use the initial conditions 

«o(0) + /3o(0) = (0ix2, 1)^, ai(0) + /?i(0) = (0ix2, 0)^. 



Then we get 



•^33 



0lx2 
1 



I 

k 



$ 



(1) 

J3 



$(1) 
^33 



fc2 



4g' 



+ 



<^;3^(o) 




-Sik^t 



(B.5) 
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From the first column of tlie equation ( IB. 21) we get 

$2x2t 






^3jt - om -±-3^- = -AeU^2x2 + 2ik{U^^2x2 - '^H^^Sj) + 4|m|2[/$2x2 + U^x'^2x2- 



$3,v - 82P<I>3 



(B.6) 



Suppose 




iUth 



Ut) , Ut) 



k 



k^ 



f---) + (^^o(t)4 



z/i(t) , z/aW 



A; 



A;2 



f--- e' 



SJfc^t 



(B.7) 

wfiere tlie coefficients C,i{t) and z^;(t), / > 0, are independent of k. 
To determine these coefficients, we substitute the above equation into 



equation (IB.6P and use the initial conditions 



?o(0) + //o(0) = (l2x2,0 



2x2,U2xi; , 



Then we get 





(B.8) 



Sik^t 



So, from the equation fIB.ll) and the asymptotic of Sj3{k) and S33(A;), 
we get the asymptotic behavior of Cj{t, A;) as A; — )■ oo. 



Cj{t,k) 



$ 



(1) 
i3 



*;? 



k^ 



+ ■■ 



(B.9) 
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